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ABSTRACT 

We present an observational evidence for the vorticity effect on the nonlinear 
evolution of the galaxy angular momentum. We first calculate the vorticity as the 
curl of the peculiar velocity field reconstructed from the 2MASS redshift survey. 
Then, measuring the alignments between the vorticity and the tidal shear fields, 
we study how the alignment strength and tendency depends on the cosmic web 
environment. It is found that in the knot and filament regions the vorticity 
vectors are anti-aligned with the directions of the maximal volume compression 
while in the void and sheet regions they are anti-aligned with the directions of 
the minimal compression. Determining the spin axes of the nearby large face-on 
and edge-on disk galaxies from the Seventh Data Release of the Sloan Digital Sky 
Survey and measuring their correlations with the vorticity vectors at the galaxy 
positions, we finally detect a clear signal of the alignments between the galaxy 
spin and the local vorticity fields. The null hypothesis that there is no alignment 
between them is rejected at the 99.999% confidence level. Our result supports 
observationally the recently proposed scenario that although the galaxy angular 
momentum originates from the initial tidal interaction in the linear regime its 
subsequent evolution is driven primarily by the cosmic vorticity field. 

Subject headings: cosmology:theory — large-scale structure of universe 



1. INTRODUCTION 

The most prevalent theory for the origin of the galaxy spins is the tidal torque theory 
which explains that a galaxy acquired its spinning motion (i.e., the angular momentum) 
via the tidal interaction with the surround i ng m atter fluctuations at its proto-galactic stage 
( IPeeblesI Il969l ; iDoroshkevichl Il970l ; White! I1984J ) . One of the key predictions of the tidal 
torque theory is the existence of the local c orrelations between the galaxy spin and the 



tidal shear fields f lCatelan fc Theunslll996al ]bl: iLee fc Penll2000l ; IPen et al.ll2000l ; iLee fc Pen 
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20011 ; IPorciani et al.l l2002al Jbl): The initial tidal torques which lasted till the turn-around 
moments caused the alignments between the spin axes of the proto-galaxies and the inter- 
mediate principal axes of the linear tidal fields. In the subsequent evolution the initially 
established spin-shear alignments would become weaker as the non-linear eff ects like galaxy 



galaxy interactions should play a role of randomiz i ng th e galaxy spin axes ([P ubinski 



2002 



20081 ). 



Catelan fc Theuns 1996b; Crittenden et al. 2001 
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Plenty of obse rvational evid ences have been found for the existence of t he intrinsic spin- 
shear correlations (lSchaferll2009l . for a comprehensive review). For instance. ILee fc Erdogdu 
(120071 ) provided a direct observational evidence for the existence of the intrinsic spin-shear 
correlations by using the tidal shear fields reconstructed from the 2MASS redshift survey 
(jErdogdu et al.l 120061 ) and the galaxy catalog complied by B . Tully f jPen et al.ll2000l ). How- 
ever, even though the main result of ILee fc Erdogdu! (120071 ) was seemingly consistent with 
the prediction of the linear tidal torque theory, they noted that the detected environmental 
dependence of the spin-shear alignments could not fit well into the tidal torque picture. The 
nonlinear effects which would destroy the initially established spin- shear alignments should 
be stronger in denser environments. Therefore, in the tidal torque picture, the galaxies 
located in denser regions were expected t o show weaker spin-shear alignments. In con- 
trast, what was found by ILee fc Erdogdu! (120071 ) was an opposite trend that the galaxies 
loca ted in denser enyiron ments exhibit stronger intrinsic spin-shear alignments (see Figure 
6 in lLee fc Erdoedul 12007b . 



A similar trend was also detected by iLed (120111 ) who measured the spin-spin alignments 
by using the nearby disk galaxies from the seventh data release of the Sloan Digital Sky 
Survey (SDSS DR7) (Ubazaiian et al.l l2009h . They showed that the galaxies which have 
ten or more neighbors within separation distance of 2 h~ 1 Mpc exhibit stronger spin-spin 
alignments. If the int rinsic spin-spin alignments were related to t he tidally induced spin 



shear alignments (e.g.JPen et al.ll2000l ; ICrittenden et al.ll2001 



2002 



Jindl2002l ; Mackev et al. 



20021 ; IPorciani et al.ll2002al Jbl). this observational finding of ILed ( 2011 ) should be interpreted 



as another evidenc e for the stronger spin-s hear alignm ents in dens er environments, consistent 
with the result of ILee fc Erdogdu! (120071 ). Although ILed (1201 if ) attempted to explain the 
observed environmental dependence of the intrinsic spin-spin alignments in the tidal torque 
picture, ascribing it to the dependence of the tidal strength on the local density, it remains 
unanswered why and how the initially established spin-shear alignments survive better in 
denser environments. 



Several literatures have also suggested that the tidal torque theory alone could not 
provide a fully satisfactory explanation for the observed dependence of the galaxy intrinsic 
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alignments on the cosmic web environment (see, e.g., 


Ha 


m et al. 


2007b: 


Jones et al. 


2010; 


Codis et al. 


2012; 


Libeskind et al. 


2012a: 


Giahi & Schafer 


2013; 


Trowland et al. 


2013) 


. The 



general consensus reached by the previous numerical works is that the observed galaxy 
intrinsic alignments are not just weakened version of the initial spin-shear correlations but 
evolved version driven by some other mechanism than the tidal shears. What has so far 
suggested as a possible candidate for the required mechanism in the literatures includes the 
large-scale cosmic flow, peculiar motions, satellite infalls, bulk flows and etc. 



Very recently, iLibeskind et al.l ( 12012b! ) proposed an interesting new scenario that the 
cosmic vorticity drives the nonlinear evolution of the galaxy angular momentum after the 
turn-around moments. In the linear regime the peculiar velocity field is curl free as it is 
described as a gradient of the perturbation potential. In the subsequent nonlinear evolution, 
however, it develops the curl mode which could affect the galaxy angular moment um. By 
analyzing the data from high-resolution N-body simulations, ILibeskind et al.l (j2012bl ) demon- 
strated that the halo angular momentum vectors are strongly aligned with the local vorticity 
vectors, indicating the presence of the dominant effect of the vorticity field on the no nlin- 
ear evolution of the galaxy angular momentum vectors. The ILibeskind et al.l (j2012bl ) also 
showed that the vorticity vectors are strongly anti-aligned with the major (minor) principal 
axes of the local tidal shear fields in the knot (void) environments. 

We note here that the puzzling environmental dependence of the galaxy spin alignments 
can be coherently explained by assuming the vorticity driven evolution of the galaxy angular 
momentu m. The goal of this Pap er is to look for an observational evidence for this new 
scenario of ILibeskind et al.l (l2012bt ). To achieve this goal, we reconstruct the cosmic vorticity 
and the galaxy spin fields from the large galaxy surveys and investigate if the vorticity-spin 
and vorticity-shear alignments found in N-body simulations exist in the real Universe. The 
upcoming three sections present the following highlights, respectively: the procedure to 
reconstruct the cosmic vorticity field and the result on the vorticity-shear alignments in 
section [2j a brief review of the algorithm to measure the galaxy spin axes with high accuracy 
and the result on the vorticity-spin alignments in section [3] a concise summary of the final 
results and discussion of the cosmological implication of our results in section HI 



VORTICITY-SHEAR ALIGNMENTS 



Erdogdu et al.l (120061 ) reconstructed the density contrast 5(x) and peculiar velocity v(x) 



fields by applyi ng the Wiener recon struction algorithm to the galaxy data from the 2MASS 
redshift survey (IHuchra et al.l 120051 ) under the assumption of a ACDM universe. The fields 
were provided on 64 3 pixels in a regular cubic space of linear size 400 /i _1 Mpc, for which the 



-4- 



position vectors x were expressed in the supergalactic coordinate system and the peculiar 
velocities v were measured in the cosmic microwave background (CMB) frame. Given the 
spatial resolution of 6.25 h~ l Mpc, the reconstr ucted fields represent ed the density and ve- 
locity fluctuations on the quasi-nonlinear scale. lErdogdu et al.l (120061 ) proved the robustness 
of their reconstruction procedure by demonstrating that the predicted density field recovers 
well the observed large scale structures of the Universe. 



Using the density field £(x ) recon structed by lErdogdu et al.l (120061 ) from the 2MASS 
redshift survey, iLee fc Erdogdul (120071 ) calculated the tidal fields as l^(x) = d^V -1 ^). 
Given that on those pixel points more distant than 100 /i~ 1 Mpc from the center the recon- 
structed density a nd velocity fields suffe red from large uncertainties (private communication 
with P.Erdogdu), lErdogdu et al.l (120061 ) selected only those 32 3 pixels whose separation dis- 
tances from the center are less than 100 h^ 1 Mpc and determined the major, intermediate and 
minor principal axes of I#(x) at each selected pixel point. The detailed descriptions of the 
re constructed density, v el ocity and tidal fields, a nd t he 2MASS red s hift su rvey can be found 



in 



Erdoedu et al.l (12009 1. ILee & Erdogdul (l2007h and lHuchra et al.l (|2005[ ). respectively. 



Now, we attempt to reconstruct the cosmic vorticity field on the selected 32 3 pixel 
points by calculating the curl of the peculiar velocity field as w = V x v. We first perform 
the Fourier transform of the peculiar velocity field to obta in its Fourier amp litudes v with 
the help of the Fast Fourier Transformation (FFT) code (IPress et al.lll992l ). The Fourier 
amplitude of the vorticity field can be written as w = k x v where k is the wave vector in 
the Fourier space, and thus the three components of w are calculated in Fourier space as 



w\ = k 2 v 3 - k 3 v 2 



w 2 = hvi - hv 3 



^3 = hv 2 - k 2 Vt . 



Finally, we perform the inverse Fourier transform of w to obtain the real-space vorticity field 
w. Figured] plots the contours of the absolute magnitude of w in the supergalactic x-y plane, 
showing the deviation of |w| from zero. Recalling the fact that the peculiar velocity field 
reconstructed from the 2MASS redshift survey was filtered on the quasi- nonlinear scale of 
~ 6.25 h~ 1 Mpc, the result shown in Figure [1] implies that even in the quasi-nonlinear regime 
the velocity field is no longer irrotational, developing the curl mode. 



As done in iLibeskind et al.l (l2012bl ). we first divide the 32 3 pixels into the knots, fila - 
ments, sheets and voids according to the signs of the shear eigenvalues (IHahn et al.ll2007al ): 
The pixels at which all three eigenvalues have positive (negative) values are marked as knots 
(voids), while the pixels at which one eigenvalue is negative (positive) and the other two 
are positive (negative) are marked as filaments (sheets). Then, we calculate the alignment 
between the vorticity vector and the principal axes of the tidal tensor at each marked pixel. 
Let {ei, e 2 , e 3 } be the major, intermediate and minor principal axes of the tidal tensor at 
each marked pixel. Calculating jj, = |w ■ e$| (for i — 1, 2, 3) and binning the values of \i in 
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the range of [0, 1], we determine the probability density distribution of fi. If there were no 
alignment, the distribution p(fi) would be uniformly unity. If p(fi) increases (decreases) with 
fi, then there should be strong alignments (anti-alignments) between the vorticity vectors 
and the principal axes of the tidal tensors. We repeat the same calculation for the knot, 
void, filament and sheet regions to separately determine p{y) for each case. 

Figure [2] plots the probability density distributions , p(/i), wit Poisson errors for the 
knot regions, showing how the vorticity vectors are aligned with the major, intermediate and 
minor principal axes of the tidal tensors in the left, middle and right panels, respectively. 
In each panel the horizontal dotted line indicates the uniform distribution of \i for the case 
of no vorticity-shear alignments. As can be seen, the vorticity vectors are strongly anti- 
aligned with the major principal axes of the tidal shear tensors, preferentially lying in the 
plane spanned by their intermediate and minor principal axes in the knot regions. In other 
words, in the highly dense knot environments the vorticity vectors are anti-aligned with the 
directions of the maximum matter compression. 

Figure [3] plots p(/i) for the void regions, revealing that in the void regions the vorticity 
vectors are anti-aligned with the minor principal axes of the tidal tensors, preferentially 
lying in the plane spanned by the major and intermediate principal axes, which is directly 
o pposite to the case of the knot regions. This result is consistent with the numerical finding 
of iLibeskind et al.l (l2012bl ) even though in their work the vorticity field was filtered on much 
smaller scale. Figures and H] plot p(fi) for the filament and sheet cases, respectively. As can 
be seen, the filament regions show no strong signal of vorticity-shear alignments, while for 
the sheet case is found a clear signal alignments (anti-alignments) with the major (minor) 
principal axes but no strong alignments with the intermediate principal axes. 



3. VORTICITY-SPIN ALIGNMENTS 

The galaxy spin axes are hard to determine in practice, even under the simplified as- 
sumption that the minor axes of the galaxies are aligned with their spin axes. In the mea- 
surements of the alignments between the galaxy spin axes and the local vorticity vectors, 
the largest uncertainty would come from the inaccurate determination of the galaxy spin 
axes. Therefore, it is important to select carefully only those galaxies whose spin orienta- 
tions can be determined with relatively high accuracy. It is often assumed that for the case 
of the late-type spiral galaxies who se shapes are close to circu lar thin discs their spin axes 



are orthogonal to the disc planes (IHaynes &: Giovanellil Il984l ). Provided that information 



on the position angles and axial ratios of the late-type spiral galaxies are available, their 
unit spin vectors can be determined up to the two-fold ambiguity in the sign of their radial 
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components (IPen et al.ll2000f ). The remaining uncertainty due to this two-fold ambiguity in 
the determination of the spin axes can be minimized by considering only those face-on or 
edge-on spirals. 

As it is essential for our investigation to determine the directions of the galaxy angular 
momentum vectors as accurately as possible, we restrict our analysis only to the nearby 
large late-type spiral galaxies viewed either face on or edge on. A sample of the nearby 



large late-type spiral galaxies was already obtained by iLed ( 120111 ) from the SDSS DR7. 
The sample contains those SP SS galaxies which have type Scd on the Hubble sequence 
( jHuertas-Company et al.ll201ll ). angular sizes larger than D c = 7.92 arcseconds in the red- 
shift range of < z < 0.02. The value of this size cut-off D c was imposed to remove the 
dwarfs which turned out to cause large uncertainties in the measurements of the spin axes 
due to their irregular shapes. Among the nearby large Scd galaxies in the sample of iLee 
( 120111 ). we make a further selection of only those ones which have axial ratios larger than 0.9 
(nearly face-on) or smaller than 0.15 (nearly edge-on) to minimize the uncertainties associ- 
ated with the two-fold ambiguity. A total of 585 nearby large late-type face-on (or edge-on) 
spirals are finally selected for our analysis. 



The unit spin vector of each selected galaxy, t = (t x , t y , i z ), is determined as ( ILee 



2011) 



: ±cos£ sin(7r/2 — 8) cos a + | sin£| sinPcos(7r/2 — 5) cos a — | sin£| cos P sin a, (2) 
y — ± cos£sin(7r/2 — 5) sina + | sin£| sinPcos(7r/2 — 5) sina + | sin£| cosPcosa, (3) 



±cos£cos(7r/2 — 5) — | sin£| sinPsin(7r/2 — 5), 



(4) 



where P is the position angle of each selected galaxy, and (5, a) are the declination and right 
ascension of each galaxy's position vector expressed in the equatorial coordinate system, and 
the plus and minus signs in front of the first terms in Equation fl2])-(jl]) represents the two-fold 
ambiguity mentioned above. Here, £ is the galaxy's inclination angle related to its axial ratio 
q and intrinsic flatness parameter p as cos 2 £ = (q 2 — p 2 )(l — p 2 ). For the galaxies of typ e 
Scd, the intrinsic flatness parameter has the value of p = 0.1 (IHaynes &: Giovanellil Il984 ) . 
Since the tidal shear and the vorticity fields from the 2MASS redshift survey have been 
reconstructed in the super-galactic coordinate systems, we find a supergalactic expression 
for the unit spin vector of each selected galaxy, s, through the coordinate transformation of 
s = Pt where P is the orthogonal matrix that transforms the equatorial to the supergalactic 
frames. 



By applying the Cloud-in-Cell interpolation (CIC) algorithm (|Hocknev &: Eastwood 



19881 ) to the tidal shear fields reconstructed by ILee fc Erdogdul (120071 ) from the 2MASS 
redshift survey, we determine the tidal shear tensors at the positions of the selected SDSS 



-7- 



galaxies and determine their principal axes, {ei, e2, e e }. Then, we calculate the cosines of 
the alignment angles, /x, between the galaxy spin axes and the principal axes of the local tidal 
tensors as /x = |s - e^|. Recall that there are two different unit spin vectors assigned to each 
selected galaxy which differ from each other by the sign of the radial components. As done in 



Lee fc Erdogdul (120071 ). we treat two spin vectors assigned to each galaxy as two independent 
realizations to end up having twice as many values of /x as the total number of the selected 
SDSS galaxies. Binning the values of /x and counting the number of those realizations, n^, 
belonging to each /x-bin, we finally determine the probability density distribution, p(/x), of 
the spin-shear alignments, calculating Poisson errors as l/(n^ — I) 1 / 2 associated with the 
determination of p(/x). 

Figure plots the probability density distributions of the cosines of the alignment angles 
between the unit spin vectors of the selected SDSS galaxies and the major, intermediate, and 
minor principal axes of the local tidal shear tensors with Poisson errors in the left, middle 
and right panels, respectively. As can be seen, the spin axes of the selected galaxies seem to 
be strongly anti-aligned with the major principal axes of the local tidal tensors, but aligned 
with the intermediate and the minor principal axes. That is, the spin axes of the selected 
SDSS galaxies tend to lie in the plane perpendicular to the local direction of maximum 
matter compression. 



Compare our result with that of iLee fc Erdogdul (120071 ) who found a significant signal 



of the correlations between the spin axes of the Tully galaxies and the intermediate principal 
axes of the local tidal tensors but no alignment signal with the other two principal axes. 
We believe that the difference resulted from the inaccurate m easurements of the spin axes 



of the Tully galaxies in the analysis of ILee fc Erdogdul (120071 ) which included not only the 



Scd galaxies but also the earlier type spirals with thick bulges without taking into account 
the two-fold ambiguity. What is newly found from our analysis is that the anti-alignments 
between the spin axes and the major principal axes are strongest while the alignments of the 
spin axes with the minor principal axes are as strong as that with the intermediate principal 
axes. Although the result is not completely against the tidal torque theory which predicts 
that the spin axes are preferentially aligned with the intermediate principal axes, we would 
like to see if there exists strong spin-vorticity alignments which may help explain better the 
detected spin-shear alignment tendency. 

Applying the CIC algorithm to the vorticity fields reconstructed in section |2j we cal- 
culate the local vorticity vectors at the positions of the selected SDSS disk galaxies. Then, 
we determine the probability distribution of the cosines of the angles between the unit spin 
vectors and the local vorticity vectors w in a similar manner, the result of which is plotted 
Figure [3 As can be seen, the probability density increases sharply and almost monotonically 



-8 - 



as fx increases, detecting a strong signal of the spin-vorticity alignments. We test the null 
hypothesis of no alignment (i.e., = 1) with the help of the x 2 -statistics and find that 
the null hypothesis is rejected at the 99.9999% confidence level. 



The result shown in Figure [7| is consistent with the numerical finding of iLibeskind et al. 



(j2012bl ). providing an observational support for the scenario that the galaxy angular mo- 
mentum evolves via its interaction w it h the l ocal vorticity in the nonlinear regime. Recalling 
that in the work of ILibeskind et al.l ( ]2012bl ) the vorticity field was filtered on the galactic 
scale of < 1 _1 h Mpc while in the current work the vorticity field is filtered on the much 
larger scale of ~ 6 /i -1 Mpc, we conclude that the vorticity effect wins over the tidal shear 
effect on the evolution of the galaxy angular momentum even in the quasi-nonlinear regime. 



4. SUMMARY AND DISCUSSION 

Utilizing the nearby lar ge face-on (or edge-on) late-type spiral galaxies from the SDSS 



DR7 (lAbazajian et al.l 120091) and the tidal she ar and vorticity fields reconstructed from the 



2MASS redshift survey (jErdogdu et al.ll2006l ). we have measured the vorticity-shear, the 



spin-shear and the spin-vorticity alignments. The reconstructed vorticity fields have spatial 
resolution of ~ 6 /i _1 Mpc, corresponding to the quasi-nonlinear regime. First, the vorticity 
vectors have been found to be strongly anti-aligned (aligned) with the major principal axes 
of the tidal shear tensors in the knot (void) regions, lying in the plane spanned by the 
other two principal axes. Second, the galaxy spin axes have turned out to be strongly anti- 
aligned with the major principal axes of the local tidal shear tensors, while aligned with the 
intermediate and minor principal axes. Finally, a clear signal of the spin-vorticity alignments 
has been detected, rejecting the null hypothesis of no spin- vorticity alignm ent at 99.9999% 



signific ance level. This result observationally stands by the new scenario of ILibeskind et al. 



( I2012bl ) that the tidally generated angular momentum of a galaxy subsequently evolves under 



the dominant effect of the vorticity field in the quasi-nonlinear and nonlinear regime. 

The spin-shear alignment tendency as well as its environmental dependence reported in 
the previous works can now be explained as follows. The cosmic flow in the nonlinear regime 
develops local vorticities which are preferentially inclined onto the plane orthogonal to the 
directions of either the maximum or the minimum volume compression depending on the web 
environment. The developed vorticities affect the galaxy angular momentum, modifying the 
spin-shear alignments from the initial tendency. Since it depends on the environment how fast 
the velocity field develops the vorticities and what directions the vorticity vectors get aligned 
with, the spin-shear alignments in the nonlinear regime come to take on the environmental 
dependence. The stronger spin-shear alignments found in denser environments should result 
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from the stronger vorticity effect there via which the nonlinear spin-shear alignments are 
established. 

An interesting cosmological implication of our result is that the nonlinear vorticity fields 
and the the vorticity-i nduced galaxy alignm ents might be usefu l as a probe of cosmology and 
gravity. As shown in iKitaura et al.l (120121 ) and mentioned in iLibeskind et al.l (j2012bl ) , the 
vorticity of cosmic fluid which equals zero in the linear regime grows in the nonlinear regime 
at third order. The more rapidly a cosmic fluid develops the third order nonlinearity, the 
stronger the spin-vorticity alignments become due to the earlier onset of the vorticity effect 
on the evolution of the galaxy angular momentum. In some modified gravity or dynamic 
dark energy models the high-order nonlinearity grows faster due to the presence of the fifth 
force than for the ACDM case. Thus, in these alternative models the faster growth of the 
nonlinearity would lead to the more rapid develo pment of the y orticity field, generating 
stronger spin-vorticity alignments. Very recently, iLi et al.l ( 120131 ) claimed that the small 
scale powers of the divergence of the peculiar velocities should be a more sensitive probe of 
modified gravity than the density power spectrum. Given our result, the small-scale powers 
of the absolute magnitude of the curl of the peculiar velocities might work powerfully as a 
complimentary probe. 
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Fig. 1. — Contours of the magnitudes of the real space vorticity fields reconstructed from 
the 2MASS redshift survey in the supergalactic x-y plane. 
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Fig. 2. — Probability density distributions of the cosines of the angles between the vorticity 
vectors and the major, intermediate and minor principal axes of the tidal shear tensors in 
the left, middle and right panels, respectively, for the knot regions where the eigenvalues of 
the tidal shear tensors are all positive. In each panel, the dotted line corresponds to the case 
of no alignment and the errors are Poissonian. 
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Fig. 3. — Same as Figure |2] but for the void regions where the shear eigenvalues are all 
negative. 
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Fig. 4. — Same as Figure [2] but for the filament regions that the largest and second to the 
largest shear eigenvalues are positive while the smallest eigenvalues are negative. 
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Fig. 5. — Same as Figure |2] but for the sheet regions that the smallest shear eigenvalues are 
negative while the other two shear eigenvalues are positive. 
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Fig. 6. — Probability density distributions of the cosine of the angles between the the galaxy 
spin vectors and the major, intermediate and minor principal axes of the tidal shear tensors 
with Poisson errors in the left, middle and right panels, respectively. 
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Fig. 7. — Probability density distribution of the cosines of the angles between the vorticity 
vectors and the spin vectors of the nearby large face-on (or edge-on) Scd galaxies selected 
from the SDSS DR7 with Poisson errors. 



